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ABSTRACT 


The horizontal sound speed in an area of complex oceano- 
graphic structure was described using cross sections obtained 
from six nonconsecutive monthly lines of STD observations at 
a 5.5 km sampling interval off Monterey, California. 

The sound speed field for each section was determined 
and visually analyzed. Cross-correlation functions of 
vertical sound speed gradients averaged over 2 m and 10 m 
increments were computed between stations. Cross-correlation 
coefficients between stations were computed for detrended 
sound speed profiles sampled at 2 m depth increments. 

Sound speed was an excellent descriptor of water mass 
features. On depth scales greater than 10 m, well defined 
sound speed field features showed horizontal extents of less 
than 11 km in some cases. On vertical scales of 2 to 10 m 
horizontal extents of less than 11 km were also evident. 
Sound speed profiles showing similarities on the scale of 2 


[wero m tended to occur at 27.5 to 38.5 km intervals. 
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ae )6UPURPOSE 

The purpose of this thesis is to study small scale 
spatial variations in the sound speed field and their time 
changes with approximately monthly time increments in the 
upper 500 m of an area associated with a complex frontal 
region. The size and character of the variability are shown 
by using isotachs to contour the sound speed field over 
each vertical section through that region. The relationship 
of the sound speed field to the water mass distribution is 
also discussed. Cross-correlation coefficients of sound 
speed profiles are utilized in order to define a sound speed 


variability length. 


B. BACKGROUND 
1. Sound Speed 

Sound speed is an increasing function of salinity, 
temperature and pressure. A number of equations relating 
sound speed to salinity (S), temperature (T) and pressure 
(P) have been developed. Wilson's October equation (Wilson, 
1960) was used to calculate sound speeds in this study. 
The importance of these variables to sound speed varies 
greatly. According to Wilson's October equation; at 7° C, 
a change of 1° C in temperature results in a 3.9 m/sec change 


in sound speed; a change of 1 ppt salinity results ina 1.3 
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m/sec change in sound speed; and a 1m change in depth 
results in a .017 m/sec change in sound speed (Kinsler and 
Pry, 1962). It is apparent, in view of the normal variations 
in the sea of T, S, and P, that temperature is by far the 
most significant variable in the upper ocean layers. The 
error resulting from the use of Wilson's October equation 
was not considered significant in the qualitative evaluation 
of the sound speed field or in the determination of the 
cross-correlation function and the cross-correlation coeffi- 
cient, since they both are measures of relative sound speed 
Wartability. 

Sound speed has two obvious uses in oceanography. 
The physical oceanographer utilizes sound speed as a property 
to define water masses and to describe water mass boundaries, 
and also to infer physical conditions (e.g., vertical advec- 
tion) by virtue of its spatial and temporal variations. Naval 
oceanographers apply directly their knowledge of sound speed 
distribution in describing the propagation of acoustic energy 
in the ocean. 

2. Frontal Zones 

A front is defined as the leading edge of a zone 
separating unlike water masses (Griffiths, 1965). A front 
is also identified in the literature as an oceanographic 
front, an oceanic front, an oceanic frontal zone, or a thermal 
front. The addition of a geographic name to any of the above 


indicates a specific feature at a specific location such as 
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the California Front, or Maltese Oceanic Frontal Zone. 
Fronts are characterized by marked horizontal and vertical 
gradients in temperature, salinity, sound speed and other 
properties. Fronts may be ill-defined, transient and 
associated with periodic vertical motions and displacements 
(Lafond and Lafond, 1966). Frontal zones usually occur at 
boundaries between current systems. They may be 200 to 400 
miles wide and consist.of a number of smaller fronts. These 
smaller features are individually and collectively referred 
Ee as fronts (Lafond and Lafond, 1971). 

Frontal zones are especially important because the 
most intense horizontal variability in oceanographic param- 
eters is found there. As a result of their importance many 
frontal zones throughout the world have been the subject of 
intense study. Two frontal zones under investigation include 
the California Front (Lafond et al., 1971) between the 
western boundary of the California Current and waters of the 
North Pacific gyre approximately 500 miles off the California 
coast; and the Maltese Oceanic Frontal Zone (Miller, 1972) 
lying East of Malta in the Mediterranean Sea. In both studies 
isothermal contours were used to identify and describe the 
frontal zone. In this discussion sound speed is used as the 
descriptive parameter. Disregarding the question of the 
relative worth of one parameter over the other as a descriptor, 
comparisons will be made freely between this study and the 


two referred to above. 
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Supporting evidence for the validity of these 
comparisons is derived from an analysis of the relative 
contributions of salinity and temperature to sound speed 


described by the following relation 


where C is sound speed, S is salinity, and T is temperature. 
The relative importance of temperature and salinity changes 


is given by comparison of the two terms within the brackets. 


They are evaluated from values of = and = given previously 
and from oe , Which can be determined from the slope of the 


ds 


temperature-salinity profiles representing this area's water 
masses as shown in figure 1. For the case where salinity 
makes the greatest contribution, indicated by the portion of 
the August, 1972 line below 450 meters, there results a 
value for = of 1.25° C/ppt. For these values the result is 


dee Oss G 4488)nedS:2 ; 


Therefore about 80% of the sound speed variability is due to 


temperature. Even in this extreme case, which is only locally 
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Significant, temperature is the dominant factor in deter- 
mining sound speed, although the contribution from salinity 
is not negligible. 
3. Oceanographic Structure in the Study Area 

The subject of study here is predominantly one 
vertical cross section running 50 miles West from a point 
(Latitude 36° 40'N, Longitude 122° O'W) located just off 
Monterey Bay. The area is shown in figures 2 and 3. This 
section transects the boundary of the colder less saline 
California Current water flowing equatorward and the warmer 
more saline water flowing poleward identified as the 
Peoeuntercurrent", “undercurrent” or "Davidson Current", 
depending on the season and the author. The structure is very 
complex in the region of strong shear between the opposing 
currents due to various dynamic mixing processes. Many 
authors including Sverdrup, Johnson and Fleming (1942), 
Reid (1963), Reid, Roden and Wyllie (1958), Wooster and Reid 
(1963), Wyllie (1966), Wooster and Jones (1970), Milnar (1972), 
Brown (1974), Wickham (1975), Blumberg (1975) and Greer (1975) 
have treated the structure of the California Current system 
and related subjects. Recent work by Wickham (1975), 
Blumberg (1975) and Greer (1975) indicates that the current 
has several branches or filaments and demonstrates variability 
both spatially and temporally on various scales. These 
authors also indicate that the transverse dimension of the 
current features is small,.especially in the poleward flow, 


on the order of tens of kilometers or less across. 
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irl. DATA COMPILATION 


A. DATA COLLECTION 

A line of 16 stations, identified by 300 series numbers, 
was established covering a total distance of 50 nm froma 
reference point at Latitude 36° 40'N, Longitude 122° 00'W 
and extending westward. This equates to a station spacing 
weer nm Or 5.5 km. On one oOccaSion, August, 1973, this 
line of stations was supplemented by two additional lines 
of stations identified by 100 series and 200 series numbers. 
Station locations are illustrated in figures 2 and 3. 

The original research plan called for stations to be 
occupied at monthly intervals over the period from August, 
1973 to August, 1974 using a Bisset-Berman continuous 
profiling salinity, temperature, and depth recorder (STD). 
Corroboration was planned to be acquired from concurrent 
Nansen casts. However, due to equipment problems and 
inclement weather, this goal could not be achieved totally. 
The goals were met for three vertical cross sections from 
August, 1973; and one each from October, 1973; November, 
1973; December, 1973; January, 1974; and August, 1974. Data 


from these sections are the basis for this study. 


B. DATA REDUCTION 
Output: from the STD is in the- form of analog traces of 


temperature and salinity versus depth. The salinity traces 


13 


a 


ae! - 
e a > ‘1a 


. - a a ion 
| 0 45 istieh-s ,eoolJiacse at de. it. 

? le solettuS fesee 8 oaixevon Bahia 
— ; : ; 
wo? ,4°OS “we sbhuvivged 27@ SRzag4 


a3 a-inge ala W3evsIeY Po ‘Bins 


a - arn 
jotted .tolsesio ate ev | 6et ELE Se 
si0be ot Wd heoveenedgqqwe sev saokear 


s© Are sfense ODL = 50 lizashe 


ra . 7 an 
wt 4 Spite + @)¢aeems. Ten pbte 
2 7 
i? 402 toec eLeves@ai ¢idaqoe Bae 


i298 otk, 61¢L ,fepgee 


a 


cy ld ot i>» loge. ,ytiaiise ; 


\dnauetq sew 10 ’ 


& 
» OF 


L wrewoH stand 


Teog elf? ,vxitmew Jem 

ite if’ 303 248 ese 
266 Sia bes ;cver 

st% tte ars 7 

cE 


136 °naoirsse 


=- 


“04 TIIY3GR 


were contaminated by erroneous fluctuations (spikes) of 
Salinity which had to be eliminated. These salinity spikes 
resulted from the inability of the conductivity cell, which 

is temperature dependent, to compensate properly for sharp 
temperature gradients. Salinity spike removal was accomplished 
by first comparing salinity traces taken as the instrument was 
lowered against traces taken as the instrument was raised. 

If the salinity spikes were present in both traces at the 

same depth in opposite directions, they were considered 
erroneous. The spikes were then smoothed using eye inter- 
polation. The error introduced by this technique was not 
considered significant due to the small salinity range and 

the small effect that salinity has on sound speed. 

The analog traces were digitalized for further analysis. 
Digitizing was done on a Calma Company Model 480 digitizer 
owned and operated by Fleet Numerical Weather Center, Monterey, 
California. The instrument converts an analog trace to a 
digital output expressed as inches of stylus travel from a 
reference point. Since stylus position is recorded every 
-01 inch, sampling interval becomes a function of scale. 

Based on the STD analog scales, .01 inch of stylus travel 
equates to .32 m of depth at the (0-300 m) scale and .75 m 
of depth at the (0-750 m) scale; a .005°C temperature incre- 
mene for a scale iof'5°C width; ‘and a:.002’ ppt) salinity 
increment for a salinity scale of 2 ppt width. The output 

is encoded on a seven track magnetic tape. Various codes are 


added as header information for identification and scaling. 
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A computer program DIGISTD, listed at the conclusion of 
this presentation, was utilized to read the seven track tape 
and convert inches of stylus travel to data values. 
Reconciliation of the data was accomplished at this point by 
the introduction of a constant correction for temperature 
(=208°C) and for salinity (.04 ppt) in the DIGISTD program 
based on a comparison of the STD data and independent concur- 
rent Nansen data. Output from the DIGISTD program was 
encoded on nine track tape, although the program allowed 
paper and card output as well. The output format gave 
temperature, salinity, sigma-t and sound speed as functions 
of depth at .3 m or .75 m increments listed by station and 
month. 

Subsequent data reduction was accomplished by generating 
an array of sound speed values as a function of depth and 
station number for each vertical cross section. Sound speed 
values were specified at whole number depth steps of 2 meters 
from 0 to 500 m or shallower at each station as required. 
Linear interpolation was utilized to assign a sound speed 
value at each designated depth. Interpolation was accomplished 
by the procedure below. Given a specified depth Dis two 
depth data points, D_y and Diy were selected such that the 
absolute value of |D-D_,| and |D-D,,| were both minimum and 
D Then, letting SV_ be the respective 


< -D < D, and SV, 


-l in dd: aE 
sound speed values of D_y and Diy , the sound speed value at 


depth D, represented by SV, was computed by the equation 
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A measure of the error in the interpolated value of SV was 
determined by calculating the mean absolute difference in 
the sound speeds at the two data points, D_, and D,, , 


averaged over one vertical cross section. 


N 
) ISv,, - SV 


MEAN SV ERROR < A 
= N l 


al 


Since the maximum mean absolute difference was .036 m/sec, 
the resulting error was considered insignificant. 

Holidays in the data were also filled using the same 
linear interpolation technique. Instead of specifying one 


depth, D, N depths dD were specified such that 


and 


where ov is the sound speed value computed at depth D,- In 
all cases, holidays were only localized and of limited 


extent. 
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When interpolated values of SV had been assigned to all 


depths at 2 m increments, mean SV gradients were calculated 
for 2m and 10 m intervals. The statistics of these 


gradients are found in the section on analysis. 
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IIIT. ANALYSIS 


A. SOUND SPEED FIELD ANALYSIS 
i Lntroduction 

The sound speed fields are illustrated in figures 4 
through 11. It has been stated previously that sound speed 
is used by physical oceanographers to identify water masses 
because, as a function of salinity, temperature, and pressure, 
it reflects changes in these variables. Higher sound speed 
tends to imply higher values of both temperature and salinity 
and vice versa. Since temperature is the dominant factor in 
determining sound speed, higher values of sound speed may 
reflect merely higher values of temperatures with lower or 
equal values of salinity, and conversely. 

The two basic water masses in the area are "southern" 
water and "northern" water and their identifying character- 
istics are higher salinity and temperature and lower salinity 
and temperature respectively. Sound speed features such as 
sound speed maximums and minimums reflect the presence of 
least mixed portions of the respective water masses. In 
considering sound speed as a water mass descriptor, sound 
speed variability due to pressure was ignored due to shallow 
depths, less than 500 meters, of the region analyzed. 

2. Major Sound Speed Field Features 
As stated previously, sound speed field features will 


now be discussed in terms of the features defined by Lafond 
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aad ouarond, (0966) ,211967) .(1967)>¢and (1971).°°oIn their 
discussions Lafond and Lafond defined characteristic features 
of the thermal field including ridges, maximums, minimums and 
frontal zones. Similar features were found and identified in 
the sound speed fields of this study. The following discus- 
sion considers the individual sound speed field features and 
their significance. 

A ridge structure is indicative of a minimum in the 
horizontal sound speed distribution and represents a colder 
less saline water mass. Examples of ridge features are shown 
at station 108 in figure 4 and station 316 in figure 9. The 
scales vary considerably between the two examples. Lafond 
and Lafond (1966) defined ridge features with horizontal 
Seales on the order of 15 nm by 45 m to 50 m in the vertical. 

Sound speed maximums and minimums are defined by 
closed isotachs and indicate extremes both horizontally and 
vertically. Maximums represent warmer more saline water 
masses of limited horizontal and vertical extent. Well 
defined sound speed maximums are illustrated at station 316 
of figure 7 and station 316 of figure 8. Less well defined 
examples are illustrated at stations 305, 314 and 316 of 
figure 10 and station 304 of figure 11. Similarly, sound 
speed minimums represent colder less saline water masses of 
limited extent. Sound speed minimums in this study are not 
as well defined as are maximums and ridge features. An 
example of a sound speed minimum occurs at station 108 in 


figure 8. 
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A frontal zone is characterized by sharp horizontal 
gradients and sloping and/or irregular isotachs marking the 
boundary between water masses. Good examples of frontal 
zones are illustrated between stations 210 and 202 in 
mraure 5; stations 313 and 311 in figure 7; at station 309 
in figure 8; and between stations 313 and 311, 308 and 305, 
and 305 and 302 in figure 10. Also of note are the apparently 
weaker frontal zones present during August, 1973 (figure 6) 
and August, 1974 (figure 11). These apparently weaker 
frontal zones are the result of more complete mixing of the 
water mass elements. However, apparently weaker frontal 
zones also result from sectioning through a frontal zone at 
some oblique angle or near an edge of the frontal zone and 
may only be a product of sampling. 

It is evident from the examples that frontal zones 
can separate any combination of sound speed minimums, sound 
speed maximums, and/or ridge features. The limbs of ridge 
features could also be considered frontal zones as they too 
identify water mass boundaries. The variability in the 
frontal zone dimensions and intensity plus the presence of 
multiple fronts in a frontal zone are also evident. Other 
researchers have described frontal zone dimensions varying 
Epomeo km (Miller, 1972) to hundreds of miles’ (Lafond et al., 
1971). 

A feature of specific Navy interest is the sharp 
vertical sound speed gradient evident in the upper 50 to 


100 m. This feature varies systematically from month to 
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month. It appears to approach the surface in August 
fe2oures 4, 5, 6, and 11). During the months of October 
(figure 7), November (figure 8), and December (figure 9), the 
feature migrates progressively downward. During December 
(figure 9), its intensity is greatly diminished and it is 
virtually absent during January (figure 10). This feature 
is seen world wide. An extensive body of research exists 
indicating that it is the result of a combination of heat 
flow across the sea surface and wind induced mixing in the 
surface layers. 

3. Comparison with Water Mass Analysis 

To further analyze the use of sound speed as a water 
mass descriptor, comparisons were made between the results of 
fares Study and the work of Lt. R. E. Blumberg (1975). 
Blumberg, using the same data, delineated water masses by 
analysis of temperature distribution relative to sigma-t 
surfaces. In his analysis Blumberg defined northern water 
as low temperature, low salinity water and southern water as 
high temperature, high salinity water. The cores of these 
water masses equate to sound speed minimums and maximums 
respectively. 

Comparisons of sound speed field features and 
Blumberg's analysis yielded close correspondence in a majority 
of cases, as illustrated in the following examples. The 
northern water centered around station 108 and the southern 
water at station 111 in figure 12 correlate well with the 


ridge feature centered at station 108 and the dip in isotachs 
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at station lll in figure 4. The northern water at station 314 
in figure 15 and station 310 in figure 16 show excellent 
agreement with the ridge feature at station 314 in figure 7 
and 310 in figure 8. The southern water at station 315 in 
figure 17 is in agreement with the dipping isotachs at 
station 315 in figure 9. The southern water at station 305 
and the northern water at station 303 in figure 18 are in 
excellent agreement with the sound speed maximum and dipping 
isotachs at station 305 and the ridge feature at station 303 
shown in figure 10. 

One point of interest that was demonstrated by a 
comparison of the two studies was the ambiguity in the origins 
of some sound speed field features. Some are derived from 
horizontal water mass variability and others arise from 
vertical motion. An illustration of this is found at 
stations 306-305 of figure 9 which indicate a sound speed 
maximum and minimum respectively. Comparison with station 
306-305 of figure 17 indicates the feature is the result of 
vertical motion in the water column and not due to the water 
mass structure. Vertical motion in figure 17 is indicated 
by the displacement of the isopycnals and isotherms in step 
with each other. 

4. Spatial and Temporal Relationships 

Inspection of the sound speed field revealed systematic 
relationships among various identifiable features. Comparison 
of the three lines of stations for August, 1973 (figure 4, 5, 


and 6) did not show the continuity of any uniquely identifiable 
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water mass feature through the three sections. This indicated 
either that the feature's spatial extent was less than the 
spacing between sections, 10 nm, or that it transited the 

area at some oblique angle. This last possibility is 
reinforced by the flow patterns shown by Wickham (1975) and 
Greer (1975), who described the region's currents using 

drogue measurements and geostrophy, respectively. 

On other vertical cross sections certain similar 
features appear on succeeding sections. A comparison of the 
300 series sections for August, 1973 and August, 1974 
(figure 6 and 11), indicated the presence of similar sound 
speed maximums, indicating a warm more saline water mass, 
at 400 to 450 meters depth in both sections. The positions 
of the features differed by 33 km between the two succeeding 
August sections. The feature in the August, 1974 section 
also appeared more well developed than its counterpart of 
the previous year. A comparison of October (figure 7), and 
November (figure 8), indicated the presence of similar well 
developed sound speed maximums representing a higher tempera- 
ture more saline water mass at station 316 in both sections. 
The feature in the November section appeared to have decreased 
definition and a greater depth by approximately 100 m than 
that in the October section. Similar ridge features 
indicating low temperature lower salinity water also appear 
in both October and November sections. Their positions vary, 


wath. one located at station 314 in.the October section and 
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in the November section. The November feature also appears 
more well developed than its October counterpart. Frontal 
zones present in both October and November sections also are 
Similar. The frontal zones in the November section are 
greatly reduced in intensity from the October frontal zones. 
The position of similar frontal zones also varied between 
sections by 11 km, being centered at station 311 in October 
and station 309 in November. A general comparison of October, 
November, December, January, the 300 series section for 
eeeuce, 1973, and August, 1974 (figures 7, 8, 9, 10, 6, and 
11) demonstrates the difficulty of predicting short-term 
changes. No obvious similarity of features exists between 
the vertical cross sections for the months of August and 
October; November and December; and December and January. 
The pair of vertical cross sections for October and November 
and for August, 1973 and August, 1974 are the only sections 
showing obvious similarity of structure. 


5. General Comments Concerning the Analysis 
of Sound Speed Features 


In concluding the discussion of the analysis of the 
sound speed fields, it is appropriate to make a few general 
comments pertaining to the overall analysis. First, visual 
scanning indicates that sound speed fields define the 
"Character" of the water mass more sharply than does the 
temperature or salinity alone. "Character" refers to the 


Shape, definition of the core and extent of the water mass. 
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Water mass elements which are small in spatial extent and 
have a definite shape and a well defined core region are 
meeustrated sat station 316 in figure 7 and at station 305 

in figure 10. Water mass structures which are large in 
extent, have less well defined cores or cores below 500 meters, 
and lack definite shape are illustrated at station 310 in 
feoure 8, Station’ 316 in “figure 9, and station 315 in 

figure 10. The extent of mixing is also evidenced by vari- 
ations in the character of the water mass structure, 
variations in the slopes of the isotachs, and the overall 
complexity of the sound speed field. The loss of the 

1487 m/sec isotach in the core region at station 316 from 
October to November (figure 7 and 8), a reduction in the core 
definition, an overall lessening of isotach slopes, and a 
reduction in the complexity of the sound speed field over the 
two months illustrate this. 

Several comments pertaining to the problems of acoustic 
propagation are also appropriate. Since the sound speed field 
ioameme controlling factor in the propagation of acoustic 
energy, its complex structure is significant. Sound speed 
minimums tend to channel acoustic energy decreasing trans- 
mission loss locally. Several poorly defined sound speed 
Minimums are illustrated at station 110 and at station 108 
in figure 4, and throughout the vertical dimension of the 
ridge feature at station 310 in figure 8. Lafond and Lafond 
(1971) found much more pronounced minimums, in terms of 


temperature, in their analysis of the California front. 
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Sound speed maximums, on the other hand, cause local 
divergence of acoustic energy. Well defined examples of 
sound speed maximums are illustrated at station 316 in 
figure 7 and station 305 in figure 10. It is apparent that 
the sound speed field in this region is very complex and 
attempts to describe it would require numerous samples of 
the vertical sound speed profile at suitable intervals in 


time and space. 


aa CROSS—-CORRELATION FUNCTIONS AND COEFFICIENTS 
i. Description of Method 

Cross-correlation functions and cross-correlation 
coefficients were computed to provide a more objective measure 
of the small scale sound speed field variability, including 
characteristic "correlation lengths". Two different methods 
of analysis were employed, one involving sound speed, the 
other its gradient. 

One method utilized the sound speed gradient averaged 
over 2 m or 10 m intervals to derive a cross-correlation 


manction (Rxy). As defined in Bendat and Piersol (1971), 


where N is the total number of depth data points; n refers 
to a specific depth data point; and x and y are sound speed 
gradients at depth n and station X and Y, respectively. The 


two Stations are separated by some horizontal distance from 
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where g; is sound speed gradient, and Ve and Va41 are sound 


speeds at depths Zs and Zend respectively, had the same effect 
as using a high pass filter. Smaller scale variations in the 
sound speed profile were emphasized and larger scale variations 
were de-emphasized. Rxy provides a numerical value repre- 
senting the small scale similarity in shape between the sound 
speed profiles at stations X and Y. 

The other method treated detrended sound speed values 
at 2 m depth increments to determine a cross-correlation 


coefficient RHOxy. As defined by Bendat and Piersol (1971) 


N 
(==) (y--¥) 
RHO my el: ae 2 a 
a ON N N 
ee —, 2 
he (xe-Rhsls sleltyszy) 
= n=l 


where N is the total number of depth data points; nis a 
specific depth data point; x and y are detrended sound speed 
values at depth n and stations X and Y respectively, separated 
by some horizontal distance from 0 to 16 station intervals; 
and x and y are the means of the detrended sound speed values 
at station X and Y respectively. Detrending was accomplished 


using the relation 
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v(n,X) = V(n,X) - yea_ .) Vy (n) 
where v(n,X) is the detrended sound speed value at a specified 
depth n and station X; V(n,X) is the sound speed value at 
depth n and station X; and the last term is the mean value 
V computed at a given depth n from the sound speed values Vv; 
at each station i=l to NSTA inclusive. Detrending also 
produced the effect of a high pass filter. Again detrending 
emphasized small scale variations and de-emphasized large 
scale variations. In this example comparing the results 
from using the sound speed gradient and the results from 
ising the detrended sound speed values would be equivalent 
to comparing the output of two high pass filters with the 
latter having alower cutoff frequency. Both methods emphasize 
small scale variations but, between the two methods, detrending 
emphasizes larger scale variations. RHOxy gives a numerical 
value between -1 and +l representing the small scale similarity 
between the shapes of the sound speed profiles at station X 
and Y. Normalization allows meaningful comparisons between 
RHOxy values determined different stations and months. 

Results are displayed utilizing two different methods. 
In one method, Rxy and RHOxy are plotted as a function of 
distance between station X and station Y. A reference station 
(Y) is designated and identified along the Y axis of each 


graph. RHOxy and Rxy are then plotted at the correlated 
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station on the X axis. Each figure contains eight graphs, 
one on top of the other. Scales, with values ranging from 


= to... 03 sect, are not indicated on the graphs 


.001 sec” 
of Rxy since their purpose is to show relative values of the 
cross-correlation function between adjacent stations. Values 
of RHOxy vary from -1 to +l as indicated. Graphs of Rxy and 
RHOxy are shown as figures 19 through 22 and figures 23 
through 28 respectively. 

Another method is also used to display the cross- 
correlation coefficients, RHOxy. Ann xn matrix of the 
cross-correlation coefficients, RHOxy, is established with 
each RHOxy being computed from a specified station X and 
station Y as indicated. Contours of equal cross-correlation 
coefficient are constructed and the resulting fields appear 
in figures 29 through 31. A maximum value of unity appears 
along the diagonal of the resulting field, with symmetry 
about the diagonal. 

2. Results of Cross-Correlation Determinations 

The cross-correlation functions derived from the sound 
speed gradients averaged over 2 meter depth increments for 
the month of October, 1973 were computed. The results are 
shown in figures 19 and 20. These cross-correlation functions 
were characterized by small correlation even between adjacent 
stations. The cross-correlation function was also computed 


from sound speed gradients averaged over 10 meter depth 


increments for October, 1973. The results are shown as 


dqe.w ddgaa le Se = re 

“Te? poatpgeet ao le? chiw xaleo’, -vaedee ahd ‘got: 

wiscwy wile #O tind gt 2h ose one 8. ea" = 
io ¢t#lev ovideley wode oF #1 Stegugg od E 
inapatbs sseweed nod tren? ao ks. i 

iS it _fetooaivat we f+ o@ I- mow3 vase 4 

7.7? bes 15 <gerowd? G1 sougptd se emote. z 

- Yisv osgese aS 


14 
iqevit of been o@fe ei Santen —tedJor 


xiogesl 3 x 3 WA .@ROET ,adasisiiisos a 
vsftellcedea St. yROGt ,ateniolitecs aciisiseke 
‘ole So/dioeqe s word bejogneo paisd iam 
. lugge lw exoeteco  .begenibak oe 
“4 

ae es2 «42 bos Batogzdenoo sae tas 
q 4 w Sa 

4 4 icy mimmanp A hf to geratt es 


id ast aie 24 lanope tS 


iJoas Sottelezaco-aaor say 
\ 

27 bopectova atnelibare 

TCL .tedessaO 40 dam 

ait O35 bag Ol ¢esupii ag 

elevyon Lieme gil Bepisetoasak 

bw oon ixiwiyves-se0u0 eT 1866 
i ie 

¢ Staelbere beoga Bavoei 


-'T°L ,aenSoteO 2oek. eftne 


#5 


figures 21 and 22. A comparison of the two cross-correlation 
functions derived from the sound speed gradients at 2 meter 
increments and 10 meter increments respectively indicated 
that better correlation existed between stations when sound 
speed gradients were averaged over 10 meter increments. This 
would be expected since the larger averaging distance would 
smooth out much of the variability due to small vertical 
motions. 

The cross-correlation coefficients derived from the 
detrended sound speed values were computed for the months of 
October, 1973, January, 1974, and August, 1974. The results 
are shown as figures 23 through 28. A comparison of these 
results and figures 7, 10, and 1l respectively indicated a 
high correlation between features in the sound speed field 
and variations in the cross-correlation coefficient. Of 
particular interest was the correlation at fronts where 
cross-correlation coefficients varied from plus values 
through zero to negative values in the distance of three 
stations, 16.5 km, or less. The algebraic difference in the 
cross-correlation coefficients over this distance exceeded 
one frequently. 

Overall the correlation of adjacent stations is small. 
In most cases the cross-correlation coefficients dropped to 
0.5 or less within a distance of two stations from the 
reference. The region of highest correlation appeared to 
be between station 303 and station 307 during October 


corgures 23, 24 and 29). 
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Correlation length, as defined here, is the distance 
between successive maximums of the cross-correlation 
coefficient. This parameter specifies the distance between 
stations having similar characteristics on the vertical 
scale of two meters or so. Correlation lengths were on the 
order of five to seven station intervals, 27.5 km to 38.5 km, 
euerng August, 1974. Correlation lengths for January, 1974 
were quite variable ranging from a distance of 5 stations, 


meoekm EO 10 stations, 55 km. 


31 


—— > 
Joatetelh off Gs , 62a Gee 


a ’ - Lor 
cotjelorey-aaase, Sie Bey 

' 

reeQtad oct 


océetele emt 


As 


a a iQ a 


a 


> a 
r" ~_ 7 
sv at? so euételsopanade deftate yi 
: i . —< - nf ie ih 

A> Oo sise aol sotte[s3365 tt 39 aisten OW. 
we 7 a7 es 

Sf of ood ©. 0°94, eiaoeetad epiteda caves od vi2© 
.Vieucst 20! appre doljelese® .6TOF «Jet 


next osiones sidsizay | #3 
tl 22 .omokseve OF ¢ 


to sorav236 « 


i 


- 


IV. CONCLUSIONS 


A. DISCUSSION 

This discussion has addressed the problems of horizontal 
scund speed variability. This was accomplished through 
subjective analysis of the sound speed field and through 
computation of cross-correlation functions and cross- 
correlation coefficients of sound speed gradients and 
detrended sound speed values respectively, for stations along 
a vertical cross section. The results, indicating a complex 
picture of sound speed variability as a function of both 
time and space, are of interest and concern to physical 
oceanographers and specifically to naval oceanographers. 

Physical oceanographers, it has been shown, can use sound 
speed as a descriptor of water masses. Its use enhances 
differences between water masses as compared to the use of 
Salinity or temperature alone where changes in salinity and 
temperature between water masses are of the same sign. For 
this reason sound speed permits high resolution in defining 
the structure of water masses. 

The naval oceanographer is interested in sound speed as 
an oceanographic variable which is the controlling factor in 
the propagation of acoustic energy. The complex sound speed 
fields in this study area result in equally complex fields 


of locally varying acoustic intensity. The nature of the 


effect of the sound speed field on the propagation of acoustic 
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energy depends on the frequency of the acoustic energy. For 
frequencies such that the wavelength is much less than the 
scale of the sound speed field feature, refraction occurs 

in a predictable manner which is not dependent on frequency. 
For frequencies whose wavelengths are on the order of the 
Size of the sound speed field feature or larger, scattering 
eccurs as a function of frequency. This would mean, for 
sound speed field features on the scale of two meters, that 


scattering becomes important near the frequency 


and, for sound speed field features on the scale of 10 meters, 


at the frequency 


where f is the frequency, C is sound speed and L is wavelength. 
It is apparent that the influence of complex sound speed 
fields, like those treated in this thesis, on acoustic 
propagation should be studied. A first step in such a study 
would be an objective description of the environment. This 
thesis provides such a description. First, for vertical 
scales of variability greater than 10 meters, it was demon- 
strated by visual analysis of the sound speed field that the 


horizontal extent of such features was, in some cases, less 
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than 11 km. Second, for scales of vertical variability 
between 2 and 10 meters, it was demonstrated by the weak 
correlation between adjacent stations that the horizontal 
extent of these features was again less than 11 km when the 
variability was integrated over the entire sound speed 
profile. The correlation length provided a measure of how 
often sound speed profiles similar on the variability scale 
of 2 to 10 meters repeated themselves. Correlation lengths 
on the order of three to five station intervals (16.5 km 


to 27.5 km) were predominant. 


B. POSSIBLE FUTURE ANALYSIS 

The application of other analysis techniques to the 
present data is suggested. One such technique is to lag the 
cross-correlation coefficient in the depth direction to 
investigate the effect of internal wave activity. Another is 
Calculation of the cross-correlation coefficients over smaller 
segments of the sound speed profile, as opposed to correlation 
over the entire depth range, to localize the centers of 
variability or homogeneity. A third variation would be to 
compute the cross-correlation coefficients over several 
different vertical sampling intervals of depth, in addition 
to the 2 and 10 meter intervals used here, to study the 
cross-correlation coefficient as a function of the variability 
scale. Future research projects might include extension of 
the present area of study to include stations further east and 
west; occupying stations at time intervals of a few days or a 


week to describe short term temporal variations; and sampling 


ae Shorter horizontal intervals. 
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Figure 1, Tyoical T-S relations in the study area (from Wickham,1975). 
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Figure 2, Area of Survey (After Wickham,1975) 
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